<?php

namespace Ycr\Pack\helpers;

/**
 * 计算距离
 * Class Distance
 * @package Ycr\Pack\helpers
 * @author ycr <Email: yuncran@126.com>
 */
class Distance
{

    /**
     * @doc计算两点地理坐标之间的距离
     * @param $longitude1 *起点经度
     * @param $latitude1  *起点纬度
     * @param $longitude2 *终点经度
     * @param $latitude2  *终点纬度
     * @param $unit       *单位 1:米 2:公里
     * @param $decimal    *精度 保留小数位数
     * @return float
     * @author ycr <Email: yuncran@126.com>
     */
    public static function getDistance($longitude1, $latitude1, $longitude2, $latitude2, $unit = 2, $decimal = 2)
    {

        $EARTH_RADIUS = 6370.996; // 地球半径系数
        $PI = 3.1415926;

        $radLat1 = $latitude1 * $PI / 180.0;
        $radLat2 = $latitude2 * $PI / 180.0;

        $radLng1 = $longitude1 * $PI / 180.0;
        $radLng2 = $longitude2 * $PI / 180.0;

        $a = $radLat1 - $radLat2;
        $b = $radLng1 - $radLng2;

        $distance = 2 * asin(sqrt(pow(sin($a / 2), 2) + cos($radLat1) * cos($radLat2) * pow(sin($b / 2), 2)));
        $distance = $distance * $EARTH_RADIUS * 1000;

        if ($unit == 2) {
            $distance = $distance / 1000;
        }

        return round($distance, $decimal);
    }


    /**
     * @doc 求两个已知经纬度之间的距离,单位为千米
     * @param $lng1 经度
     * @param $lat1 纬度
     * @param $lng2
     * @param $lat2
     * @return float  距离，单位千米
     * @author ycr <Email: yuncran@126.com>
     */
    public function getDistanceTwo($lng1, $lat1, $lng2, $lat2)//根据经纬度计算距离
    {
        //将角度转为弧度
        $radLat1 = deg2rad($lat1);
        $radLat2 = deg2rad($lat2);
        $radLng1 = deg2rad($lng1);
        $radLng2 = deg2rad($lng2);
        $a = $radLat1 - $radLat2;//两纬度之差,纬度<90
        $b = $radLng1 - $radLng2;//两经度之差纬度<180
        $s = 2 * asin(sqrt(pow(sin($a / 2), 2) + cos($radLat1) * cos($radLat2) * pow(sin($b / 2), 2))) * 6378.137;
        return $s;
    }

    /**
     * select * from sb_user_geo order by
     * ACOS(SIN(('.$city_lat.' * 3.1415) / 180 ) *
     * SIN((lat * 3.1415) / 180 ) + COS(('.$city_lat.' * 3.1415) / 180 ) * COS((lat * 3.1415) / 180 ) *
     * COS(('.$city_lng.' * 3.1415) / 180 - (lng * 3.1415) / 180 ) ) *
     * 6380  asc  limit 30
     */

}